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NOTATION 


A  Panel  area 

[  A  |  Matrix  of  circulation 

A,(  Influence  coefficients 

|C|  A  vector  defined  in  Eq  (3.6) 

Cn  Local  normal-force  coefficient  for  a  fixed  axial  location  on  the  body 
Cp  Pressure  coefficient 

d.  Velocity  induced  normal  to  control  point  i 

ds  Area  element 

dT  Volume  element  associated  with  the  vorticity  wT(q) 

G  Loop  circulations 

|  G  |  Loop  circulations 

i  Dummy  variable 

k  Numbers  of  iteration 

I,  Vectors  connecting  consecutive  nodes  of  the  panel 

m  Number  of  panel  edges 

N  Number  of  panels  on  the  body 

n  Unit  normal 

P  Local  pressure 

Poo  Freestream  dynamic  pressure 

q  Position  vector 

r  Radius  of  the  body 

r  Position  vector 

SDM  Source-distribution  method 

s  A  scalar 

V  Surface  velocity 
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^  Velocity  vector 


^oo  Undisturbed  freestream  velocity 

Vco  Undisturbed  freestream  velocity 

VLM  Vortex  -  lattice  method 

w  Vorticity  vector 

ws  Surface  vorticity  distribution 

u  T  Vorticity  vector 

T  Circulation 

(  Representative  thickness  of  each  region  of  nonzero  vorticity 
0  Circumferential  angle  measured  from  the  windward  meridian 

q  Fluid  density 

r  Volume  of  the  vortex  field 


ABSTRACT 


A  nonlinear  three-dimensional'  vortex-lattice  method  was  developed  which 
treats  the  steady  separated  flow  over  prolate  bodies  with  open  separation 
moving  through  an  inviscid  incompressible  fluid.  The  strength  and  position  of 
the  body  wake  were  found  as  part  of  the  solution.  Specifically,  flows  with 
smooth-surface  separation  were  considered  as  opposed  to  flows  with  sharp-' 
edge  separation  treated  with  the  vortex-lattice  concept  in  the  past.  To 
demonstrate  the  technique,  results  for  the  flow  over  an  inclined  ogive-cylinder 
are  presented.  - 

In  the  case  of  attached  flow,  comparisons  are  presented  of  the  results 
from  the  vortex-lattice  method  using  optimal  and  average  control  point  loca¬ 
tions  with  the  results  of  the  source-distribution  method  and  with  experimental 
data.  The  same  panel  arrangement  was  used  in-the  calculations  for  both 
methods.  The  results  demonstrated  that  the  results  of  the  present  method  are 
somewhat  more  sensitive  to  panel  arrangement  than  are  those  of  the  source- 
distribution  method.  Also,  the  effect  of  control  point  location  varies 
dramatically  as  the  incidence  of  the  body  is  changed. 

In  the  case  of  separated  flow,  results  of  the  vortex-lattice  method  were 
compared  with  experimental  data  and  with  the  results  of  a  typical  two- 
dimensional  analogy.  The  results  demonstrated  that' the  present  method  agrees 
most  favorably  with  the  experimental  data  windward  of  a  separation  line. 


ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  performed  under  the  Naval  Sea  Systems  Command  Special  Focus 
Program  on  Separated  Flow  Field  on  Submarines  during  High-Speed  Turns,  and  was  authorized  by  the 
General  Hydrodynamics  Research  Program  administered  by  the  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center  (DTNSRDC).  Funding  was  provided  under  Program  Element  61153N,  Task  Area 
SR02VM0B.  and  Work  Unit  1542-102. 


ACKNOWLEDGMENTS 

I  would  like  to  acknowledge  the  encouragement  and  friendship  of  my  colleagues  from  the  Ship  Perfor¬ 
mance  Department.  Hydromechanics  Branch.  Their  contributions  proved  invaluable.  Specifically,  I  thank  my 
supervisors.  Dr.  Thomas  Huang  and  Justin  McCarthy  for  their  support  and  technical  expertise.  In  addition.  I 
am  particularly  grateful  to  Professor  Mook  for  his  unsurpassable  assistance  provided  during  my  graduate  and 
undergraduate  years  at  Virginia  Polytechnic  Institute. 


1 .0  INTRODUCTION  AND  LITERATURE  REVIEW 


The  problem  of  potential  flow  separating  from  a  smooth  surface  is  fundamentally  important  in  fluid 
dynamics  In  fact,  one  encounters  smooth-surface  separation  in  the  potential  flow  over  almost  any  body 
executing  some  maneuver  of  interest.  Examples  of  problems  with  smooth-surface  flow  are  the  flows  about: 

•  Missiles  at  angles  of  attack 

•  Wings  at  high  angles  of  attack  or  yaw 

•  Aircraft  fusilages 

•  Ship  stems 

Ke.eav.  .-“numerical  schemes  based  on  the  vortex-lattice  concept  have  had  considerable  success  in  treating 
the  lull)  three  dimensional  separated  potential  flow  over  low  aspect-ratio  thin  wings  at  high  incidence  and 
the  attached  potential  flow  over  bodies  of  more  general  shape.  So  far,  separation  has  been  constrained  to 
take  place  along  a  sharp  edge.  However,  in  principle,  the  vortex-lattice  method  is  not  limited  in  application 
so  long  as  convection  dominates  vorticity  diffusion.  / 

We  developed  a  nonlinear  three-dimensional  vortex-lattice  method  that  treats  the  steady  incompressible 
iiniseid  flow  over  bodies  in  cases  where  the  separation  takes  place  along  a  smooth  surface.  We  assumed  that 
the  flow  does  not  reattach  and  vortex  breakdown  does  not  take  place  near  the  body.  The  location  of  the 
separation  lines  on  the  body  were  assumed  known.  The  strength  and  shape  of  the  body  wake  were  found  as 
part  of  the  solution. 

I'o  demonstrate  the  feasibility  of  the  technique,  we  present  results  for  the  flow  over  an  inclined  body 
with  a  tangent ogive  shaped  nose  and  a  cylindrical  afterbody,  hereafter  referred  to  as  an  “ogive-cylinder.” 

I  he  separation  line  used  in  the  calculations  begins  aft  of  the  ogive-cylinder  junction  and  has  a  fixed  cir- 
cumtcrential  location  along  the  length  of  the  body. 

In  the  remainder  of  this  section  we  discuss  the  flow  topology  and  review  similar  efforts  to  solve  the 
problem 


I  I  I  <  >[’<>[  (Hill  \l.  DISCUSSIONS 

f'lie  lust  fundamental  work  concerned  with  the  topological  notions  used  here  was  that  of  Poincare  (1928) 
but  ilu-.  work  followed  the  more  recent  (1980)  development  of  Peake  and  Toback.  We  considered  the  flow 
about  a  smooth  slender  body  inclined  slightly  to  a  uniform  oncoming  freestream.  One  of  the  streamlines  in 
the  oncoming  flow  attaches  itself  to  the  body  near  the  nose  at  the  forward  stagnation  point.  The  attachment 
point  aits  as  a  source  ot  skin-friction  lines  which  emerge  and  envelop  the  body.  All  of  the  skin-friction  lines 
end  at  the  rear  stagnation  point,  which  acts  as  a  sink  for  the  lines.  From  this  rear  stagnation  point  a  single 
streamline  enters  the  external  flow  and  extends  downstream.  The  forward  stagnation  point  is  referred  to  as  a 
nodal  singular  point  ot  attachment,  and  the  rear  stagnation  point  as  a  nodal  singular  point  of  separation. 

An  attached  flow.  then,  was  defined  as  a  flow  containing  two  and  only  two  nodal  points  in  the  pattern 
ot  skin  friction  lines.  These  are  the  nodal  points  of  attachment  and  separation.  A  three-dimensional  flow  with 
more  than  two  smeul.u  nodal  points  in  the  skin-friction  line  pattern  was  defined  as  a  separated  flow'. 

flic  simplest  three  dimensional  separated  flow  contains  three  nodal  points,  two  being  of  the  same  type. 
Suppose  that  there  are  two  nodal  points  of  attachment  and  one  nodal  point  of  separation.  The  skin-friction 
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lines  emerging  from  the  two  attachment  points  must  be  prevented  from  crossing.  This  requires  another 
singular  point  between  them,  commonly  known  as  a  saddle  point.  One  skin-friction  line  from  each  of  the 
attachment  points  enters  the  saddle  point  and  two  skin-friction  lines  emerge.  These  emerging  skin-friction 
lines  act  as  separatrices  in  the  skin-friction  line  pattern.  The  other  skin-friction  lines  tend  to  converge  along 
either  side  of  the  separatrices.  The  separatrices  are  the  separation  lines  on  the  body  and  they  always  emerge 
from  a  saddle  point.  The  skin-friction  line  patterns  of  attached  flows  never  contain  any  saddle  points, 
whereas  the  skin-friction  line  patterns  of  separated  flows  always  contain  at  least  one  saddle  point.  The  con¬ 
vergence  of  skin-friction  lines  along  either  side  of  a  particular  skin-friction  line  does  not  necessarily  con¬ 
stitute  separation:  a  skin-friction  line  emanating  from  a  nodal  point  cannot  be  a  separation  line. 

In  this  case  there  are  three  nodal  points  and  one  saddle  point,  which  makes  the  difference  between  the 
number  of  nodal  and  saddle  points  the  required  two.  As  a  general  rule  the  number  of  nodal  points  in  the 
skin-friction  line  pattern  must  exceed  the  number  of  saddle  points  by  two  (Davey,  1961). 

A  surface  extends  from  the  body  into  the  outer  flow.  It  originates  along  a  separation  line  and  prevents 
the  boundary  layer  on  cither  side  of  it  from  merging.  The  surface  rolls  up  as  it  extends  downstream  and  is 
the  body  wake.  Vorticity  generated  within  the  boundary  layer  on  the  body  is  shed  into  the  wake  along  a 
separation  line. 

Sears  (1948)  was  among  the  first  to  introduce  the  concept  of  “limiting  streamlines”  while  discussing  the 
laminar  boundary  layer  on  an  infinite  yawed  cylinder.  If  we  let  the  height  of  a  streamline  in  the  boundary 
layer  on  the  body  approach  zero,  the  normal  velocity  tends  to  vanish  faster  than  the  tangential  velocity. 

Thus,  the  resultant  streamline  or  limiting  streamline  is  parallel  to  the  body  and  its  projection  onto  the  body 
surface  is  coincident  with  a  skin-friction  line. 

Limiting  streamlines  have  behavior  similar  to  that  of  skin-friction  lines;  that  is.  they  emerge  from  nodal 
points  of  attachment,  disappear  into  nodal  points  of  separation,  and  converge  along  certain  limiting 
streamlines  that  are  parallel  to  separation  lines.  Sears  included  a  sketch  indicating  this  behavior  near  a 
separation  line  but  did  not  refer  to  it  as  such.  Eichelbrenner  (1954)  and  Eichelbrenner  and  Oudart  (1955) 
proposed  that  a  three-dimensional  separation  line  was  the  envelope  of  converging  adjacent  limiting 
streamlines. 


FREE  VORTEX  LINES 


OPEN  SEPARATION  CLOSED  SEPARATION 

figure  I  I  Schematic  Views  of  Open  and  Closed  Separation  Patterns 
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Each  curved  vortex  line  in  the  wake  is  replaced  with  a  number  of  short  straight  connected  vortex 
segments.  The  first  segment  in  each  wake  line  attaches  to  the  bound-vortex  lattice  and  the  rest  extend 
downstream  a  Unite  distance.  At  the  end  of  the  last  segment  for  each  line  we  place  a  straight  semi-infinite 
vortex  line  which  extends  downstream  parallel  to  the  undisturbed  freestream. 

The  separation  line  condition  demands  that  the  bound  vorticity  coincident  with  a  separation  line  vanish. 

In  view  of  the  role  that  the  longitudinal  bound-vortex  segments  play  in  modeling  the  bound  vorticity,  we  set 
the  circulation  around  each  of  the  longitudinal  vortex  segments  along  a  separation  line  to  zero.  Spatial  con¬ 
servation  of  circulation  determines  the  strength  of  the  free-vortex  segments.  The  wake  lines  are  included  in 
ihe  calculation  of  the  influence  coefficients,  which  correspond  to  the  panels  that  border  a  separation  line. 

This  makes  the  influence  coefficient  matrix  a  function  of  the  wake  shape. 

The  problem  remains  to  determine  the  orientation  of  the  first  wake  segment  in  each  line.  J.H.B.  Smith 
(1978)  has  shown  that  a  free  sheet  leaves  the  body  tangentially  as  the  Reynolds  number  approaches  infinity. 
We  are  unable  to  allow  the  sheet  to  do  this,  however,  due  to  numerical  difficulties  with  the  singular  nature 
of  vortex  lines  I'o  circumvent  this  problem  at  present  we  simply  place  the  first  segment  in  each  wake  line 
perpendicular  to  the  body  and  make  it  short. 

The  Kelvin  Helmliot/  theory  of  vorticity  requires  that  free  vorticity  travel  downstream  with  the  local  par¬ 
ticle  velocity.  Thus,  each  wake  vortex  line  is  actually  a  streakline.  If  we  view  the  steady  state  solution  as  a 
typical  time  step  m  the  solution  of  the  more  general  unsteady  problem  with  steady  flow  conditions,  each  of 
the  finite  wake  vortex  segments  represents  an  Euler  step  in  the  time  integration  of  the  path  followed  by  a 
fluid  particle.  We  align  each  finite  wake  vortex  segment  with,  and  make  its  length  proportional  to,  the  local 
particle  velocity.  This  provides  tor  a  force-free  wake  The  proportionality  constant  is  the  time  increment, 
here  chosen  to  be  unity  (he  velocity  at  the  upstream  end  of  a  wake  segment  is  used  for  directing  and  sizing 
the  segment.  Kelly  ( 1977)  investigated  the  use  of  other  points  along  the  segment  and  found  that  the  fewest 
iterations  were  required  using  the  upstream  end  and  that  the  predicted  loads  were  best  as  well.  Schroder 
(1978)  used  the  average  ot  the  velocities  calculated  at  the  upstream  and  downstream  ends,  but  this  practice 
almost  doubles  the  computer  time  necessary  for  solution  and  the  results  were  not  significantly  different  from 
those  obtained  by  using  the  upstream  end 

.V4  THE  SEPARATED  1  1  <>W  SOLUTION  PROCEDURE 

The  task  is  to  find  the  position  ot  each  finite  vortex  segment  in  the  wake  and  to  find  the  set  oi  loop  cir¬ 
culations  that  simultaneously  enforce  the  no-penetration  condition  on  the  body  and  render  the  wav.  force- 
tree  The  follow  mg  iteration  procedure  is  used: 

1  Present  the  wake  position 

2  Calculate  the  influence  coefficient  matrix. 

7  Solve  tor  the  loop  circulations  which  enforce  the  no-penetration  condition. 

4  Position  the  wake  vortex  segments  so  that  they  are  force-free. 

V  Repeat  steps  2  to  4  until  convergence  is  achieved. 

In  step  2.  we  set  the  wake  position  either  to  be  flat  or  to  be  obtained  from  a  previous  solution.  Step  2 
involves  changing  only  those  columns  in  the  influence  coefficient  matrix  which  correspond  to  panels  border¬ 
ing  a  separtion  line  Steps  t  arid  4  require  more  elaborate  calculations  and  are  discussed  in  Sections  3.b  and 
3.7 


The  influence  coefficient  A,j  is  equal  to  the  normal  velocity  induced  at  control  point  i  due  to  a  unit  loop 
circulation  around  panel  j.  The  solution  of  Equation  (3.3)  is  addressed  in  Section  3.6. 

In  the  case  of  attached  flow,  once  the  no-penetration  condition  is  satisfied  we  can  calculate  the  loads  on 
the  body,  as  described  below  in  the  section  on  calculation  of  loads. 

3.3  THE  FREE- VORTEX  LATTICE 

The  free-vortex  lattices  replace  the  continuous  free-vortex  sheets  in  much  the  same  way  that  a  bound- 
vortex  lattice  replaces  the  bound-vortex  sheet.  Each  free-vo.  ;.'x  sheet  joints  the  bound-vortex  sheet  along  a 
prescribed  separation  line. 

First  arrange  the  bound-vortex  lattice  such  that  each  separation  line  coincides  with  one  of  the 
longitudinal  vortex  lines  on  the  body.  Furthermore,  the  circumferential  vortex  lines  are  placed  so  that  each 
separation  line  on  the  body  begins  and  ends  on  a  node.  Such  an  arrangement  is  possible  only  when  the 
separation  line's  circumferential  variation  is  small  The  arrangement  is  shown  schematically  in  Figure  3-4. 


Figure  3-4  -  Schematic  View  of  the  Attachment  of  a  Free- 
and  Bound-Vortex  Lattice 


A  number  of  semi-infinite  nonintersecting  curved  vortex  lines  are  placed  on  each  of  the  free-vortex 
sheets  Each  wake  vortex  line  emanates  from  a  node  on  the  bound-vortex  lattice  on  a  separation  line  and  ex¬ 
tends  downstream  parallel  to  a  streamline.  These  wake  vortex  lines  divide  each  free-vortex  sheet  into  stream- 
wise  strips  or  "ribbons''  of  surface  vorticity.  The  vorticity  distributed  on  each  of  these  wake  ribbons  is 
subsequently  concentrated  into  the  vortex  lines  that  border  the  wake  in  a  fashion  similar  to  the  way  that  the 
bound  vorticity  is  concentrated  into  the  bound-vortex  segments. 
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Spatial  conservation  of  vorticity  (or  alternatively,  circulation)  demands  that  the  sum  of  the  branch  cir¬ 
culations  at  any  node  in  a  vortex  lattice  vanish.  A  convenient  way  to  satisfy  this  requirement  is  to  define  the 
branch  circulations  in  terms  of  loop  circulations.  This  is  analogous  to  the  familiar  loop  and  branch  currents 
used  in  the  analysis  of  electrical  networks.  Denote  loop  circulations  with  the  letter  G.  Figure  3-3  shows  that 


the  branch  circulation  1~  for  the  vortex  segment  bordering  panels  labeled  i  and  i+  1  is  given  by: 

r  =  Gi+I-  Gj  (3.2) 

We  solve  for  the  loop  circulations  as  primary  variables,  then  immediately  write  the  branch  circulations  as 
algebraic  combinations  of  the  loop  circulations.  The  Biot-Savart  law  is  used  to  calculate  the  velocity  field 
due  to  the  individual  vortex  segment  with  its  associated  branch  circulation.  (See  section  on  velocity  calcula¬ 
tions  below.)  To  find  the  velocity  at  a  field  point  we  add  the  freestream  velocity  to  the  combined  velocity 
due  to  all  the  vortex  segments  in  the  flow. 


3.2  THE  ATTACHED  FLOW  SOLUTION  PROCEDURE 


To  enforce  the  no-penetration  condition,  a  set  of  simultaneous  linear  equations  is  written  for  the  loop  cir¬ 
culations,  or: 

N 

£  AjjGj  =  -^*n,  fori  =  1.2 . N  (3.3) 

J=1 


where  i 

A,i 

Voo 

n. 

and  N 


denotes  the  panel  number. 

are  the  influence  coefficients. 

is  the  undisturbed  freestream  velocity. 

is  the  unit  normal  for  panei  i. 

is  the  number  of  panels  on  the  body. 
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Figure  3-2  -  Representation  of  Surface  Vorticity  by  Bound-Vortex  Segments 


m 


r  ,i,a 


A  ~  (A  +  A) 
i  =  1  i 


(3.1) 


where  i  is  the  side  of  the  panel. 

1 1  denote  the  vectors  connecting  consecutive  nodes  of  the  panel, 

T,  denotes  the  branch  circulations  of  the  corresponding  vortex  segments, 

A  denotes  the  panel  area. 

A,  denotes  the  area  of  the  adjacent  panel  on  the  i-th  side, 
m  is  the  number  of  panel  edges  (either  three  or  four). 

Branch  circulations  used  in  the  pressure  calculations  that  are  calculated  with  this  area  weighting  are 
referred  to  as  area-weighted  circulations.  Branch  circulations  calculated  using  equal  contributions  from  each 
of  the  neighbonng  panels  all  referred  to  as  equally-weighted  circulations. 

As  demonstrated  in  the  next  section,  weighting  the  branch  circulations  helps  alleviate  some  of  the  pro¬ 
blems  normally  encountered  with  irregularly  spaced  lattices. 

In  this  formulation  the  vorticity  strength  is  constant  across  each  panel  and  is  discontinous  along  panel 
edges.  A  higher-order  formulation  would  allow  the  vorticity  strength  to  vary  across  each  panel  surface  and 
be  continuous  to  some  degree  along  the  panel  edges. 
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Short  curved  vortex  segments  that  connect  the  various  nodes  on  the  body  act  as  edges  of  the  vortex 
panels.  We  replace  these  curved  segments  with  straight  segments  that  connect  the  same  nodes.  The  panels 
now  have  straight  sides.  These  straight  vortex  segments  make  up  the  bound-vortex  lattice. 

The  three-sided  vortex  panels  are  flat.  The  four-sided  vortex  panels  are  nonplanar,  in  general,  because 
all  four  comers  need  not  lie  in  the  same  plane.  In  practical  applications,  however,  the  four-sided  panels  are 
almost  Oat.  Thus,  the  panel  area  and  normal  direction  for  a  nonplanar  panel  are  calculated  as  though  the 
panel  were  Oat.  The  normal  direction  of  a  panel  is  parallel  to  the  cross  product  of  the  panel’s  diagonals  and 
its  area  is  equal  to  one-half  the  magnitude  of  the  same  cross  product. 

The  no-pentration  condition  is  satisfied  on  the  body  using  a  collocation  procedure.  That  is,  the  no¬ 
penetration  condition  is  enforced  at  a  finite  number  of  discrete  points  on  the  body  and  not  elsewhere.  (In¬ 
deed,  large  unrealistic  normal  velocities  are  induced  between  collocation  points.)  These  collocation  points  are 
called  “control  points"  and  one  is  placed  on  each  panel. 

The  panel's  control  point  is  placed  on  a  plane  that  contains  the  centroid  of  the  comers  of  the  panel  and 
is  perpendicular  to  the  panel  s  normal.  This  plane  is  referred  to  as  the  panel  plane.  If  the  panel  is  nonplanar 
then  the  panel's  centroid  usually  will  not  lie  on  the  panel  surface  nor  will  the  panel  plane  contain  any  of  the 
panel’s  comers. 

In  their  distributed  source  method,  Hess  and  Smith  (1962)  defined  a  panel  plane  for  each  panel  in  the 
same  manner  used  here  and  then  project  the  panel’s  comers  onto  the  panel  plane  along  a  direction  parallel  to 
the  panel's  normal.  Thus,  with  their  method  the  comers  of  adjoining  panels  are  not  necessarily  coincident  as 
they  are  here. 

Two  choices  for  the  placement  of  the  control  point  on  a  panel  plane  were  explored.  The  first  choice,  us¬ 
ed  by  most  authors,  is  to  place  the  control  point  at  the  panel's  centroid.  This  set  of  control  points  is  called 
“average"  control  points.  A  second  choice,  investigated  by  Kelly  (1977)  and  used  by  Asfar.  Mook,  and 
Nayfch  (1978),  is  to  place  the  control  point  at  the  point  on  the  panel  plane  of  minimum  normal  velocity  in¬ 
duced  by  a  unit  loop  circulation  (defined  below)  around  the  panel.  These  control  points  are  called  "optimal” 
control  points,  and  their  use  can  have  a  profound  effect  on  the  results  obtained,  as  shown  in  the  next 
section. 

The  circulation  around  a  vortex  segment,  called  a  “branch"  circulation  is  denoted  with  the  symbol  f.  Its 
sign  and  associated  direction  is  taken  according  to  the  right-hand  rule. 

Each  vortex  segment  replaces  the  surface  vorticity  component  parallel  to  it  and  in  an  area  immediately 
surrounding  it.  Thus,  longitudinal  segments  replace  the  longitudinal  components  of  vorticity  on  the  body  sur¬ 
face  and  likewise  for  the  circumferential  direction.  Figure  3-2  illustrates  this  concept.  The  longitudinal  com¬ 
ponent  of  vorticity  in  shaded  area  “A"  is  replaced  by  the  vortex  segment  “a"  while  the  circumferential 
component  of  vorticity  in  shaded  area  “B”  is  replaced  by  the  vortex  segment  “b.”  If  the  two  panels  border¬ 
ing  a  vortex  segment  have  equal  areas  then  the  segment  receives  one-half  of  its  strength  from  each  of  those 
panels.  That  is.  the  vortex  segment's  strength  originates  equally  from  each  of  the  panels  the  segment 
borders.  However,  if  the  two  panels  have  areas  which  are  widely  different  then  more  of  the  segment's 
strength  originates  from  the  panel  with  the  larger  surface  area.  In  an  attempt  to  account  for  this  discrepancy, 
the  branch  circulations  are  weighted  with  the  appropriate  panel  areas  when  approximating  the  vorticity  vec¬ 
tor,  w,  associated  with  a  vortex  panel  or: 


13 


3  0  NUMERICAL  MODELING 


In  this  section,  a  numerical  model  and  solution  procedure  are  presented  to  accompany  the  hydrodynamic 
model  presented  in  Section  2.0.  The  model  is  based  on  the  vortex-lattice  concept  of  replacing  continuous 
vortex  sheets  w  ith  arrangements  of  vortex  segments.  The  placement  of  the  vortex  segments  on  the  vortex 
sheets  is  guided  by  the  physics  of  the  How.  We  replace  the  bound-  and  free-vortex  sheets  with  lattcies  of 
short  straight  vortex  segments.  The  position  of  the  segments  in  the  free-vortex  sheet  and  the  strength  of  the 
segments  in  both  the  bound-  and  free-vortex  sheets  are  then  found  by  a  nonlinear  iteration  procedure. 

First  the  modeling  of  the  bound-vortex  sheet  is  described  and  the  procedure  for  solving  for  the  attached 
tlow  over  the  body  explained.  Next,  the  modeling  of  the  free-vortex  sheets  is  described  and  the  iteration  pro¬ 
cedure  outlined  for  solving  for  the  separated  tlow  over  the  body.  We  then  discuss  velocity  calculations,  the 
solution  of  the  arising  set  of  simultaneous  linear  equations,  the  iteration  process  to  render  the  wake  force- 
free.  and  finally  the  calculation  of  the  loads  on  the  body.  Many  of  the  concepts  introduced  in  this  section  are 
illustrated  further  in  the  next  section  with  numerical  examples. 

3  1  THE  BOUND-VORTEX  LATTICE 

We  begin  by  replacing  the  continuous  bound-vortex  sheet  with  circumferential  and  longitudinal  families 
of  curved  vortex  lines  which  lie  on  the  body's  surface.  The  longitudinal  vortex  lines  start  at  the  nose  of  the 
body  and  continue  back  to  the  body's  tail.  The  circumferential  vortex  lines  form  a  series  of  axially  spaced 
vortex  rings  which  enclose  the  body  These  two  families  of  vortex  lines  intersect  at  various  nodes  to  form 
rows  and  columns  of  vortex  panels  or  elements  which  cover  the  entire  body.  The  vortex  lines  are  placed 
such  that  the  resulting  arrangement  of  vortex  panels  is  regular  and  uniform.  The  panels  that  meet  to  form  the 
nose  and  tail  are  three-sided  and  the  remaining  panels  are  four-sided.  Such  an  arrangement  is  shown 
schematically  in  Figure  3-1. 


Figure  3-1  -  Schematic  View  ol  a  Bound-Vortex  Lattice 
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2.3  SEPARATION  LINE  CONDITIONS 


We  must  specify  conditions  along  the  separation  lines  as  well.  Since  we  are  posing  our  problem  in  terms 
of  a  vorticity  distribution,  we  seek  a  condition  on  the  strength  of  vorticity  along  a  separation  line.  J.H.B 
Smith  ( IdXO)  showed  that  the  inviscid  streamlines  just  leeward  and  windward  of  a  separation  line  are  parallel 
and  inclined  to  the  separation  line,  respectively.  The  windward  streamlines  then  leave  the  body  surface  and 
vorticity  generated  in  the  boundary  layer  is  swept  into  the  outer  flow  along  them.  It  follows  that  nonzero 
bound  vorticity  along  a  separation  line  would  generate  a  crossflow  relative  to  the  separation  line  on  the  sur¬ 
face  of  the  body.  Since  the  streamlines  just  leeward  of  a  separation  line  are  parallel  to  the  separation  line, 
we  set  the  bound-vortex  strength  coincident  with  and  parallel  to  the  zero  separation  line.  Spatial  conservation 
of  vorticity  demands  that  this  vorticity  be  swept  into  the  outer  flow  where  it  travels  downstream  to  form  the 
body  wake. 

2.4  SUMMARY 

In  this  section  we  have  discussed  the  approximation  of  steady  separated  flow  over  a  body  by  a  con¬ 
tinuous  vortex-sheet  model.  Solving  for  the  shape  of  the  free-vortex  sheets  and  the  vorticity  distribution 
everywhere  remains  subject  to  the  following  conditions: 

1.  There  is  no  penetration  through  the  bound-vortex  sheet. 

2.  Vorticity  is  divergenceless  or  alternatively  there  is  spatial  conservation  of  circulation. 

3  Free  vorticity  is  convected  downstream  at  the  local  particle  velocity. 

4.  Bound  vorticity  coincident  with  a  separation  line  vanishes.  The  requirement  of  spatial  conservation 
of  circulation  then  dictates  the  strength  of  the  wake. 

In  the  next  section  we  present  a  numerical  model  for  the  hydrodynamic  model  presented  here. 


where  q  is  the  position  vector  locating  the  volume  element  dT  with  associated  vorticity  wT(q),  the  integration 
is  performed  over  the  volume  of  the  vortex  field,  t,  and  the  curl  is  taken  with  respect  to  the  components  of  r. 

The  regions  of  nonzero  vorticity  are  thin  everywhere.  This  leads  us  to  let  the  representative  thickness  of 
each  of  these  regions,  i,  tend  toward  zero  and  simultaneously  to  let  the  strength  of  the  vorticity  tend  to 
infinity  in  such  a  manner  that  the  far-field  influence  remains  the  same.  While  taking  the  limit,  the  product  of 
the  thickness  and  the  strength  of  the  vortex  sheet  remains  constant.  Thus,  the  boundary  layer  on  the  body 
becomes  a  bound-vortex  sheet  and  the  thin  wake  regions  each  become  free-vortex  sheets.  We  define  the  sur¬ 
face  vorticity  distribution.  ws,  on  these  vortex  sheets  as: 

ws(q)  =  lim  £  wT(q)  (2.2) 

£  ♦  0 
wT»  00 


We  now  integrate  liquation  (2.1)  over  the  thickness  and  obtain: 


*  ♦  l 
V(  r)  =  — 
4rr 


(2.3) 


where  the  position  vector  q  now  locates  the  area  element  ds  and  the  integration  is  performed  over  the 
combined  surfaces  of  the  bound-  and  free-vortex  sheets.  The  problem  remains  to  specify  the  shape  of  the 
free-vortex  sheet  and  the  vorticity  distribution  on  the  bound-  and  free-vortex  sheets. 


2.2  KINEMATIC  FLOW  CONDITIONS 

There  are  a  number  of  conditions  that  the  solution  must  satisfy.  First,  there  must  be  no  penetration 
through  the  body  surface.  We  write  that: 


^  •  n  =  -  'O'oo  •  n  (2.4) 

where  ’0r«>  is  the  undisturbed  freestream  velocity.  ^  is  the  disturbance  velocity,  and  n  is  the  unite  normal  to 
the  body  surface. 

Second,  vorticity  is  a  divergenceless  field  by  construction  or: 

div(w)  =  div  (curl  ^)  =  0  .  (2.5) 

Furthermore,  according  to  the  Kelvin-FIelmhotz  theory  of  vorticity,  the  substantial  derivative  of  circula¬ 
tion.  T.  on  a  free-vortex  sheet  is  zero  or: 

DT  =  0  (2.6) 

Dt 

Consequently,  free  vorticity  is  conceded  downstream  with  the  local  particle  velocity.  Thus,  a  steady  free- 
vortex  sheet  is  a  stream  surface. 
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separation  to  take  place  along  sharp  edges.  We  compare  the  results  of  the  vortex-lattice  and  source- 
distribution  methods.  In  the  calculations,  the  same  lattice  configuration  was  used  for  both  methods.  To  our 
knowledge,  this  is  the  only  published  comparison  of  the  two  methods  where  calculations  were  performed 
using  precisely  the  same  lattice  layout.  We  also  show  that  the  effect  of  control  point  placement  varies 
significantly  with  incidence.  This  discovery  sheds  new  light  on  some  modifications  to  the  vortex-lattice 
method.  These  modifications  were  tested  only  at  zero  incidence  for  original  publication. 

Section  2.0  is  a  mathematical  statement  of  the  problem  and  development  of  a  hydrodynamic  model  based 
on  continuous  sheets  of  vorticity. 

Section  3.0  develops  a  numerical  model,  based  on  the  vortex-lattice  concept,  to  accompany  the 
hydrodynamic  model. 

Section  4.0  compares  the  results  of  the  method  with  those  of  other  methods  and  experiments  and 
investigates  several  numerical  phenomena  as  well. 

Section  5.0  contains  conclusions  about  the  work  and  recommendations  for  further  areas  of  study. 

2.0  HYDRODYNAMIC  MODELING 

Here  we  describe  the  hydrodynamic  model  for  which  a  numerical  model  was  developed.  Consider  the 
steady  inviscid  incompressible  flow  over  a  body  at  an  angle  of  attack.  The  fluid  moves  with  a  uniform 
freestream  velocity  far  from  the  body  and  its  extent  is  infinite.  The  angle  of  attack  is  sufficient  to  develop 
large  scale  separation  on  the  lee  side  of  the  afterbody  but  not  so  large  as  to  precipitate  unsteady  vortex  shed¬ 
ding.  Any  asymmetry  in  the  flow  is  due  to  geometry  rather  than  to  flow  instabilities.  Reattachment  and 
secondary  separation  effects  are  neglected. 

2.1  BOUND  -  AND  FREE- VORTEX  SHEETS 

When  the  flow  separates,  vorticity  generated  within  the  boundary  layer  on  the  body  is  shed  into  the 
outer  flow  along  separation  lines.  This  shed  vorticity  forms  the  body  wake.  The  resulting  flow  is  fully  three- 
dimensional  and  nonlinear.  The  wake  shape  and  strength  are  important  in  determining  the  flow  field 
surrounding  the  body  and  ultimately  the  loads. 

If  the  Reynolds  number  is  high,  the  wake  is  thin  when  compared  with  the  overall  dimensions  of  the 
body.  Consequently,  the  vorticity  in  the  flow  is  concentrated  into  the  thin  boundary  layer  on  the  body  and 
the  thin  vortical  regions  comprising  the  wake.  The  flow  outside  these  regions  :s  essentially  irrotational .  We 
pose  the  problem  in  terms  of  the  strength  of  the  vorticity  distribution  and  the  shape  of  the  body  wake. 

We  can  express  the  velocity  field,  V,  at  a  field  point  located  by  the  position  vector,' r,  in  terms  of  the 
vorticity  distribution,  w,.  for  any  incompressible  flow  as  (see  Karamcheti,  1966,  Chapter  9): 
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1.3.2  Separated  Flows 

The  three-dimensional  flow  over  thin  wings  with  sharp-edge  separation  has  been  treated  with  nonlinear 
vortex-lattice  methods  by  a  number  of  authors  over  a  wide  variety  of  planforms.  Steady  flows  were  con¬ 
sidered  by  Mook  and  Maddox  (1974);  Kandil  (1974);  Kandil,  Mook,  and  Nayfeh  (1976);  as  well  as  Zorea 
and  Rom  (1978).  Kelly  (1977)  presented  a  valuable  parametric  study  of  the  steady  nonlinear  vortex-lattice 
method  for  thin  wings. 

Kandil,  Chu,  and  Yates  (1980)  describe  a  ‘Hybrid’  vortex-lattice  method  whereby  a  continuous-vorticity 
model  is  used  for  near-field  velocity  calculations  and  a  discrete- vortex  model  (such  as  the  one  used  here)  is 
used  for  far-field  calculations.  The  method’s  feasibility  was  investigated  and  some  preliminary  results 
presented  but  significant  work  remained  to  be  done.  Yen,  Mook,  and  Nayfeh  (1981)  developed  a  continuous- 
vorticity  panel  method  enforcing  continuity  of  vorticity  strength  across  panel  edges.  Their  results  were 
encouraging. 

Nikolitsch  (1978)  treated  wing-body  combinations  by  combining  Wardlaw’s  version  of  Bryson’s  vortex  line 
model  with  Gersten’s  (1961)  nonlinear  lifting-surface  theory  for  wings.  Wing-body  combinations  were  also 
treated  by  Uchiyama.  Mikkilineni,  and  Wu  (1978)  using  distributed  sources  on  the  surface  of  the  body  and  a 
nonlinear  vortex-lattice  methixl  to  model  the  wings.  They  attached  a  straight  semi-infinite  vortex  line  on  either 
side  of  the  nose  of  the  body.  The  attachment  points  and  strength  of  each  of  the  vortex  lines  were  empirically 
estimated.  Sheffield  and  Deffenhaugh  ( 1980)  used  the  same  approach  to  treat  simple  bodies  at  incidence  but 
solved  for  the  shape  of  the  vortex  lines  as  well. 

Atta  and  Nayfeh  (1978)  solved  for  the  flow  over  wing-body  combinations  using  a  vortex-lattice  method 
for  both  the  body  and  the  wings.  They  allov  ed  for  separation  off  the  sharp  edges  of  the  wing  only. 

Maskew,  Rao,  and  Dvorak  ((980)  and  Maskew  ((981)  used  a  combination  doublet-lattice  and  source- 
distribution  method  enforcing  a  Dirichlet  boundary  condition  to  solve  for  the  flow  over  thick  wings  where 
flow  separation  was  not  confined  to  the  edges  of  the  wing  and  thus  treated  a  problem  similar  to  the  one  we 
consider  here. 

Fiddes  (1980)  solved  for  separation  angle  on  a  cone  at  incidence  by  interacting  F.T.  Smith’s  (1978)  ex¬ 
tension  of  the  two-dimensional  tnple-deck  boundary-layer  theory  of  Sychev  (1972)  with  a  sophisticated 
vortex-sheet  model.  Fiddes  represented  the  vortex  sheets  as  a  series  of  circular  arcs  which  meet  tangentially 
at  their  end  points.  Their  application  was  somewhat  limited  as  Fiddes  employed  slender-body  and  conical 
(low  assumptions  and  his  boundary -layer  method  was  limited  to  laminar  (lows. 

Rehbach  (1977)  presented  an  entirely  different  approach  to  represent  three-dimensional  vortex  Hows 
adopting  a  Lagrangian  view  of  the  circulation  shed  from  a  separation  line.  He  used  a  vortex  lattice  on  the 
wing  .surface  but  integrated  the  vorticity  transport  equation  in  the  wake.  There  was  some  difficulty  in 
matching  the  two  presentations  of  vorticity  along  the  separation  line,  which  he  described *as  “delicate ” 

1.4  SUMMARY 

This  section  introduced  the  topic  of  discussion  and  briefly  reviewed  the  recent  literature  pertaining  to  the 
problem.  More  complete  general  reviews  of  three-dimensional  separated  Hows  are  given  by  J.H.B.  Smith 
(1975.  1977.  i980).  Peake  and  Tobak  (1980).  and  Tobak  and  Peake  (1981a,  1981b). 

In  a  departure  from  earlier  work,  we  present  a  vortex-lattice  method  that  treats  the  flow  over  bodies  with 
separation  along  a  smooth  surface.  Farlier  vortex-lattice  methods  either  treated  attached  flows  or  constrained 
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Typically,  a  Neumann  boundary  condition  (i.e.,  no  penetration  through  the  body  surface)  is  enforced  in 
a  collocation  manner.  The  collocation  points  are  referred  to  as  “null”  or  “control”  points.  A  set  of  linear 
algebraic  equations  is  written  for  the  strengths  of  a  finite  number  of  discrete  singularities  on  or  near  the  body 
surface.  Upon  solution,  the  fluid  velocity  (or  velocity  potential)  is  evaluated  and  the  pressure  coefficient 
determined  at  each  control  point. 

One  of  the  earliest  integral-equation  methods  is  that  of  Hess  and  Smith  (1962,  1964,  1966).  Here,  a 
distribution  of  fiat  quadrilateral  source  panels  of  constant  strength  are  placed  near  the  body  surface  in  the 
manner  described  by  Kellogg  (1929).  The  strength  of  each  source  panel  is  determined,  and  the  velocity 
potential  is  calculated  at  each  control  point.  Numerical  differentiation  along  the  body  surface  yields  the  fluid 
velocity  at  each  control  point.  Bernoulli’s  equation  is  then  used  to  calculate  the  pressure  coefficient.  Land- 
weber  and  Macugno  (1969)  modified  the  method  in  the  treatment  of  the  singularity  of  the  kernel  in  the 
numerical  procedure.  Webster  (1975)  developed  a  method  using  triangular  source  panels  of  linearly  varying 
strength.  In  his  model,  the  source  strength  was  continuous  across  panel  edges.  In  an  effort  to  reduce  surface 
velocity  irregularities,  he  submerged  the  singularity  surface  beneath  the  body  surface  as  was  done  previously 
by  Pien  (1964).  The  method  calculated  relatively  smooth  velocity  variations  on  the  body  surface  and  allowed 
for  fairly  accurate  calculation  of  streamlines. 

Chang  and  Pien  (1975)  used  a  distribution  of  doublets  on  the  body  surface  in  a  method  similar  to  the 
aforementioned  source-distribution  methods.  Johnson  and  Rubbert  (1975)  combined  sources  and  doublets  in  a 
higher-order  method  which  enforced  either  a  Neumann  or  Dirichlet  boundary  condition  on  the  velocity  poten¬ 
tial  in  a  collocation  manner.  The  source  strength  for  each  panel  was  prescribed  to  cancel  the  component  of 
the  freestream  velocity  normal  to  the  panel.  They  found  that  using  a  Dirichlet  boundary  condition  resulted  in 
remarkable  insensitivity  of  the  results  to  panel  arrangement. 

Asfar  (1978)  treated  the  flow  over  bodies  of  general  shape  with  a  surface  distribution  of  vortex  segments 
(a  vortex-lattice  method).  In  the  vortex-lattice  method,  the  fluid  velocity  at  the  control  points  is  calculated 
directly  from  the  singularity  strengths  rather  than  differentiating  the  velocity  potential.  Asfar,  Mook,  and 
Nayfeh  (1978)  showed  that  combining  sources  of  known  strengths  with  the  vortex-lattice  method  resulted  in 
more  accurate  pressure  distributions  over  bodies  of  simple  shape.  In  this  case,  the  source  strengths  were 
adjusted  to  cancel  a  portion  of  the  normal  component  of  the  freestream  velocity  for  each  panel. 

In  all  the  methods  mentioned  thus  far,  the  strength  of  an  auxiliary  singularity  distribution  is  first 
obtained  by  construction  of  a  large  influence  coefficient  matrix.  Then,  the  velocity  (or  velocity  potential)  is 
calculated  on  the  body  by  solving  a  set  of  simultaneous  linear  equations.  Another  approach  is  to  solve  for 
the  singularity  distribution  on  the  body  in  an  iterative  fashion.  The  key  advantage  in  this  alternative  method 
of  solution  is  that  far  less  computer  storage  is  required  and  often  the  iterative  procedure  can  be  shown  to 
converge  rapidly  for  problems  of  interest.  Among  the  first  to  apply  the  iterative  approach  to  fluid  flow  pro¬ 
blems  were  Chow.  Hou  and  Landweber  (1976).  Chow.  Lee.  and  Owen  (1982)  used  an  iterative  approach  to 
solve  for  the  source  distribution  on  a  body  while  Noblesse  (1983)  and  Noblesse  and  Triantafyllou  (1983) 
developed  a  method  to  solve  iteratively  for  the  velocity  potential  on  the  body.  Among  the  attractive  features 
of  Noblesse's  technique  is  that  especially  rapid  convergence  was  proved  for  the  case  of  longitudinal  transla¬ 
tion  of  a  slender  or  thin  body.  It  also  can  be  extended  to  treat  frec-surface  flow  problems. 
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with  the  vortex  sweepback  principle  led  to  satisfactory  predictions  of  the  vortex  strength.  In  the  sweepback 
principle,  the  wake  is  considered  to  be  part  of  a  steady  infinite  vortex  street.  The  vortex  street  velocity  is 
then  cancelled  by  the  component  of  the  freestream  normal  to  the  vortex  lines.  In  the  application  of  the  prin¬ 
ciple,  it  is  important  that  the  cross-flow  component  of  the  freestream  can  be  considered  incompressible. 

Kubin  (1973),  Fidler  and  Bateman  (1975),  and  Kao  (1975)  use  this  approach  with  minor  variations  to 
account  for  the  nose  vorticies  which  do  not  follow  the  pattern  of  a  Karman  vortex  street.  The  Thompson  and 
Morrison  correlations  provided  the  required  empirical  data. 

While  these  two-dimensional  approaches  are  by  no  means  completely  adequate  to  model  the  complicated 
picture  of  a  general  three-dimensional  separated  flow,  it  is  intriguing  that  they  can  provide  us  with  more  than 
just  qualitative  details  of  the  flow.  Several  engineering  problems  have  been  solved  using  the  cross-flow  plane 
analogy.  It  is  perhaps  important  to  note  that  the  topologies  of  two-dimensional  unsteady  and  three- 
dimensional  steady  flows  appear  analogous  (Tobak  and  Peake,  1981a,  1981b).  Stricter  correspondence 
regarding  the  difference  between  the  two  problems  certainly  leaves  areas  of  some  doubt.  For  example,  the 
growth  of  an  unsteady  two-dimensional  vortex  differs  in  important  ways  from  the  development  of  a  steady 
three-dimensional  vortex  structure.  In  three  dimensions,  the  fluid  entering  the  core  of  the  vortex  can  be 
discharged  axially  whereas  in  two  dimensions  such  an  escape  is  not  possible.  Thus,  a  two-dimensional  vortex 
must  expand  to  accommodate  fluid  entering  the  core  region.  Only  in  one  case  are  the  two  kinds  of  vortex 
flows  formally  identical:  a  steady,  three-dimensional  conical  vortex  and  an  unsteady  two-dimensional  vortex 
that  grows  linearly  with  time.  See  Peake  and  Tobak  (1980)  for  a  discussion  of  the  issue  of  two-dimensional 
and  three-dimensional  vortex  correspondence. 

Another  important  issue  dealing  with  two-dimensional  analogs  is  the  lack  of  upstream  influence.  The 
upstream  influence  in  the  flow  is  more  pronounced  when  the  body’s  radius  varies  along  the  length  of  the 
body.  One  attempt  to  account  for  the  expanding  radius  of  the  body  as  one  moves  axially  along  the  nose  is  to 
assume  that  the  cylinder  being  modeled  in  the  impulsive  flow  analogy  is  expanding.  The  expansion  rate  is 
related  to  the  number  of  cross-flow  planes  used  in  the  calculation  and  to  the  speed  and  incidence  of  the 
body . 

1.3  THREE-DIMENSIONAL  APPROACHES 
1.3.1  Attached  Flows 

The  fundamental  problem  of  the  three-dimensional  attached  flow  over  an  arbitrary  body  has  been  treated 
by  a  number  of  authors.  The  most  popular  and  successful  of  these  methods,  have  been  integral-equation 
methods.  These  methods  express  the  solution  of  the  problem  in  terms  of  an  integral  equation  which  is  then 
treated  numerically.  They  differ  from  one  another  mainly  in  the  form  of  the  kernel  and  the  numerical  pro¬ 
cedure  used  to  obtain  a  solution.  Integral-equation  methods  effectively  reduce  the  three-dimensional  problem 
to  a  two-dimensional  one.  since  the  velocity  (or  velocity  potential)  need  be  calculated  only  on  the  body  and 
wake  surfaces,  not  in  the  complete  external  flow.  However,  only  with  the  advent  of  high-speed  computers 
have  these  methods  proved  practical. 

The  use  of  an  integral-equation  method  can  be  broken  up  into  two  main  parts.  First,  the  form  of  the 
integral-equation  form  must  be  chosen.  This  is  equivalent  to  choosing  the  type  of  singularity  (e.g.,  sources, 
doublets,  vorticies,  etc.)  Second,  the  solution  procedure  and  boundary -condition  form  (Neumann  or 
Dirichlet)  must  be  specified. 


vortex  is  introduced  in  the  outer  flow  at  each  separation  point.  The  vorticies  in  the  wake  subsequently  travel 
downsteam  under  their  mutual  influence. 

The  discrete-vortex  method  is  plagued  with  a  number  of  numerical  problems,  however.  Birkhoff  and 
Fisher  (1959)  and  Hama  and  Burke  (1960)  repeatc  .  PvOsenhead’s  calculations  and  found  that  increasing  the 
number  of  vorticies  tracked  in  the  wake  often  led  to  their  chaotic  motion.  Fink  and  Soh  (1974)  developed  a 
rediscretization  procedure  whereby  the  vorticity  in  the  wake  is  redistributed  among  the  vorticies  at  each  time 
step.  However.  Baker  (1980)  claimed  that  their  test  case  was  inadequate  and  their  error  estimate  suspect. 
Baker  used  their  technique  to  study  the  wake  behind  a  ring  wing  and  found  that  the  discrete-vortex  method 
still  predicted  physically  unrealistic  results  as  the  number  of  vorticies  was  increased  when  using  Fink  and 
Soh's  modification.  Baker  could  not  arrive  at  a  definitive  reason  for  the  breakdown,  although  it  appears  that 
he  did  not  reduce  the  temporal  step  size  in  a  consistent  manner  as  he  increased  the  number  vorticies  tracked 
in  the  wake.  In  any  case,  such  a  rediscretization  process  in  three  dimensions  would  probably  be 
unmanageable. 

There  are  a  number  of  empirical  factors  in  the  discrete-vortex  formulation  as  well.  For  example,  the 
predicted  strength  of  the  vorticies  shed  into  the  wake  must  be  reduced  significantly  for  realistic  results. 
Schoaff  (1978)  and  Sarpkaya  and  Schoaff  (1979)  decreased  the  strength  of  the  vorticies  gradually  as  they 
travelled  downstream  in  an  effort  to  simulate  the  dissipation  and  diffusion  of  vorticity.  The  shape  and  size  of 
their  dissipation  function  was  guided  by  numerical  experimentation.  Significantly,  they  solved  for  the  separa¬ 
tion  angle  numerically  and  found  it  to  agree  favorably  with  experiments. 

The  cross-flow  analogy,  in  combination  with  the  discrete-vortex  method,  has  been  used  to  calculate  the 
flow  over  bodies  at  incidence  with  varying  degrees  of  success.  Unfortunately,  the  approach  requires  em¬ 
pirical  guidance  and  considerable  artwork  on  the  part  of  the  user.  For  examples  of  the  discrete-vortex  method 
and  its  application,  see  Sarpkaya  (1966),  Gerrard  (1967),  Sarpkaya  (1968),  and  Laird  (1971). 

A  similar  approach  using  the  cross-flow  plane  analogy  uses  a  lumped  vortex  cross-flow  model  (Bryson, 
1959).  Here,  a  point  vortex  is  joined  to  the  body  with  a  feeding  sheet  of  negligible  strength.  The  system  is 
then  adjusted  such  that  it  is  globally  force-free  with  the  force  on  the  feeding  sheet  balancing  the  force  on  the 
point  vortex.  The  force  per  unit  of  length  on  the  body  is  equal  to  the  time  rate  of  change  of  the  momentum 
in  the  cross-flow  plane.  Schindel  (1965,  1969)  used  conformal  mapping  to  extend  this  method  to  bodies  with 
elliptic  cross-sections  while  Wardlaw  (1973)  calculated  side  forces  on  bodies  combining  asymmetric  vorticies 
and  slender-body  theory. 

In  a  modification  of  Bryson’s  model,  called  the  multi-vortex  model,  a  number  of  free  point  vortices  were 
allowed  to  roll  up  and  form  a  concentrated  vortex  in  each  cross-flow  plane.  Angelucci  (1971)  and  Marshall 
and  Deffenbaugh  (1974)  treated  bodies  with  symmetric  vortex  shedding  in  this  way,  while  Wardlaw  (1975) 
extended  the  method  to  treat  asymmetric  vortex  shedding.  Angelucci  calculated  the  force  distribution  using 
the  vortex  impulse  theorem  while  Wardlaw  used  a  momentum  balance  similar  to  Bryson.  Marshall  and  Def¬ 
fenbaugh  calculated  and  then  numerically  integrated  the  circumferential  pressure  distribution  to  find  the 
normal-force  coefficient  on  the  body.  They  also  solved  for  the  boundary-layer  flow  numerically  using  two- 
dimensional  criteria  to  determine  the  location  of  the  separation  lines  on  the  body. 

Thompson  and  Morrison  (1971)  found  that  analyzing  the  wake  behind  an  inclined  cone-cylinder  by 
means  of  the  impulsive  flow  analogy  led  to  reasonable  predictions  of  the  Strouhal  number,  but  the  vortex 
strength  was  over-estimated.  However,  they  found  that  the  Karman  (1930)  vortex  street  theory  combined 
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Maskell  (1955)  suggested  that  the  limiting  streamlines  on  opposite  sides  of  a  separation  line  merged  and 
left  the  surface  as  single  streamlines  on  the  “surface  of  separation”  (wake).  He  categorized  simple  separated 
flows  into  two  classes:  a  "closed”  or  bubble-type  and  an  “open”  or  shear-layer  type.  These  are  shown 
schematically  in  Figure  1-1.  If  the  body  is  symmetric,  the  saddle  point  is  located  on  the  plane  of  symmetry 
lor  closed  separation  and  off  the  plane  of  symmetry  for  open  separation.  In  the  case  of  open  separation,  fluid 
coming  from  upstream  of  the  body  wets  both  sides  of  the  wake  and  the  pressure  is  continuous  across  it. 

Legendre  (1956,  1965)  and  Oswatitsch  (1957)  independently  showed  that  by  viewing  the  pattern  of 
limiting  streamlines  as  a  continuous  vector  field,  the  behavior  of  the  merging  of  streamlines  called  for  in 
Masked's  discussion  could  not  be  accommodated.  Reproducing  their  results  and  in  a  departure  from  his 
earlier  work.  Lighthill  (1963)  chose  to  work  with  skin-friction  and  surface-vortex  lines.  Both  sets  of  lines 
cover  the  body  completely,  intersect  at  right  angles,  and  are  unique  everywhere.  Moreover,  the  separation 
line  itself  is  a  skin- friction  line.  Since  the  skin-friction  lines  are  uniquely  defined,  their  pattern  can  be  viewed 
as  a  continuous  vector  field.  He  showed  that  the  number  and  type  of  singular  points  on  the  surface  obey 
topological  rules.  However,  Legendre  (1977)  noted  that  as  far  as  was  then  currently  known,  all  experi¬ 
mentally  determined  skin-friction  line  patterns  could  be  interpreted  within  the  framework  of  his  earlier  work 
and  that  a  more  complicated  singular-behavior  theory  was  not  necessary  to  explain  observed  patterns. 

The  rules  governing  the  number  and  type  of  singular  points  can  be  extended  above  the  surface  on  planes 
of  symmetry  (Hunt,  Abell,  Peterka,  and  Woo,  1978),  on  projection  of  conical  flows  (J.H.B.  Smith,  1969) 
and  on  cross-flow  planes  (Perry  and  Fairlie,  1974). 

We  now  turn  to  a  discussion  of  other  attempts  to  treat  the  problem  of  interest. 


1.2  TWO-DIMENSIONAL  APPROACHES 

A  large  number  of  methods  treat  the  separated  flow  over  a  prolate  body  by  means  of  a  two-dimensional 
analogy.  Typically,  the  three-dimensional  steady  problem  is  related  to  the  solution  of  a  two-dimensional 
time-dependent  problem.  The  two-dimensional  unsteady  problem  is  solved  and  then  an  approximate  solution 
of  the  three-dimensional  problem  inferred. 

Many  two-dimensional  methods  employ  the  “cross-flow  plane”  or  “impulsive  flow"  analogy,  first  sug¬ 
gested  by  Allen  anil  Perkins  (1951)  These  methods  consider  the  impulsively  started  flow  over  a  cylinder  to 
he  analogous  to  the  flow  in  a  plane  perpendicular  to  the  body’s  longitudinal  axis.  Time  in  the  former  pro¬ 
blem  is  related  to  the  axial  distance  of  the  cross-flow  plane  in  the  latter  problem.  Hence,  the  task  remains  to 
find  the  two-dimensional  time  dependent  flow  over  a  cylinder  experiencing  impulsively  started  flow. 

The  most  rigorous  formulations  using  this  approach  numerically  solve  the  Navier- Stokes  equations  for 
the  impulsive  flow  over  the  cylinder.  Thoman  and  Szewczyk  (1969)  and  Jordan  and  Fropim  (1972)  treated 
the  flow  over  an  inclined  cylinder  in  this  manner  while  Walitt  and  Trulio  (1971)  considered  the  compressible 
flow  over  an  inclined  body  of  revolution.  Unfortunately,  the  computational  requirements  for  this  approach 
can  be  prohibitive  due  to  the  Navicr-Stokes  integration. 

Less  rigorous  and  more  practical  formulations  employing  the  cross-flow  plane  analogy  use  the  discrete- 
vortex  method  to  obtain  the  impulsive  flow  solution.  Originally  developed  by  Rosenhead  (1932).  the 
discrete-vortex  method  is  based  on  replacing  the  Karman  vortex  streets  behind  the  moving  cylinder  with 
point  vorticiex.  Time  is  discretized  and  marched  in  a  pedestrian  fashion.  At  each  time  step,  a  new  point 
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1.5  VELOCITY  CALCULATIONS 


The  velocity  field  due  to  a  vortex  segment  is  calculated  using  the  Biot-Savart  law  (see  Karamcheti, 

1966).  For  a  single  straight  vortex  segment  of  strength  T.  the  velocity  at  field  point  [  given  by  the  follow¬ 
ing  computationally  advantageous  form  (see  Figure  3-5): 


Biot-Savart  Law 


VP  =  - 
4  n 


hi  *  h  , 

|  hi  *  n| 

L  rajj 

(3.4) 


The  velocity  due  to  a  vortex  segment  tends  to  infinity  as  the  field  point  approaches  the  segment.  To 
avoid  this  physically  inconsistent  behavior,  we  insert  a  viscous  core  around  the  segment  where  the  induced 
velocity  is  zero. 


3.6  THE  SOLUTION  OF  SIMULTANEOUS  LINEAR  EQUATIONS 

An  efficient  solution  of  the  arising  set  of  simultaneous  linear  equations  is  important  in  solving  for  the 
flow  over  the  body.  This  is  especially  true  in  the  case  of  separated  flow  since  we  must  do  it  once  for  every 
iteration  between  the  wake  position  and  loop  circulations. 

In  general,  iterative  methods  are  particularly  well  suited  for  solving  linear  systems  of  equations  when  a 
good  initial  guess  is  available  or  when  the  number  of  equations  prohibits  its  complete  storage  in  high-speed 
memory.  For  an  iterative  method  we  use  the  modified  Gauss-Seidel  method  (Johnson  and  Riess,  1977)  or  the 
modified  Gauss-Seidel  method  with  successive  overrelaxation,  sometimes  called  the  SOR  method  (Goult. 
Hoskins,  Milner,  and  Pratt.  1974). 
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In  contrast,  direct  methods  are  particularly  well  suited  for  solving  linear  systems  of  equations  when  a  good 
initial  guess  is  not  available  or  when  the  number  of  equations  allows  for  complete  storage  in  high-speed  mem¬ 
ory.  For  a  direct  method,  we  use  the  Gauss-elimination  method  with  partial  pivoting  (Johnson  and  Riess,  1977). 

The  number  of  operations  (multiplications  and  divisions)  for  the  Gauss-Seidel  method  is  0(N2)  per  itera¬ 
tion,  while  the  Gauss-elimination  method  requires  0(N3/3)  operations  for  complete  solution.  Of  course,  with 
the  direct  method,  considerable  row  exchanges  may  be  necessary  if  pivoting  is  used.  Row  exchanges  do  not 
require  significant  additional  computation  if  row  indexing  is  used,  however. 

A  number  of  factors  determine  the  structure  and  conditioning  of  [A]  and  thus  govern  the  choice  of  solu¬ 
tion  scheme.  By  inspection,  we  can  show  that  if  the  body  is  closed,  for  any  row  i  in  [A]: 

N 

T  Ay  =  0  (3.5) 

j  =  l 

Thus,  every  column  of  (A]  is  a  linear  combination  of  the  remaining  columns  and  the  matrix  is  singular. 
This  is  because  the  loop  circulations  are  the  unknowns  in  the  problem  while  the  branch  circulations  deter¬ 
mine  the  velocity  field.  Since  each  of  the  branch  circulations  is  defined  as  the  difference  between  two  of  the 
loop  circulations,  the  set  of  loop  circulations  which  induce  a  particular  velocity  field  are  unique  only  within 
an  additive  constant  or  functionally: 

V([G])  =  V([G]  +  s[C])  (3.6) 

where  s  is  a  scalar  and  [Cj  is  a  vector  of  ones  and  is  of  the  same  length  as  [G]. 

If  the  body  is  not  closed  then  we  can  show  that  for  any  row  i  in  [A]: 

N 

^  Ajj  =  dj  (3.7) 

J»1 

where  dj  is  the  velocity  induced  normal  to  control  point  i  due  to  a  unit  branch  circulation  around  each  of 
the  vortex  segments  which  border  an  opening. 

For  those  rows  which  correspond  to  panels  far  away  from  any  opening,  dj  can  be  O(10-4)  but  it  grows 
rapidly  as  we  move  towards  an  opening.  Even  though  the  matrix  is  no  longer  singular,  partial  pivoting  was 
necessary  for  the  Gauss-elimination  method  to  obtain  a  solution. 

While  using  the  modified  Gauss-Seidel  method,  we  found  that  after  several  iterations,  successive  iterates 
of  the  loop  circulations  in  both  the  singular  and  nonsingular  cases  could  be  approximated  by: 

|G]k  +  !  =  [G]k  +  sk[C]  (3.8) 

where  the  superscript  indicates  the  iterate  number  and  sk  is  a  positive  scalar.  By  monitoring  the  convergence 
of  a  few  of  the  branch  circulations  during  the  Gauss-Seidel  iteration,  we  were  always  able  to  find  a  solution. 
In  most  cases,  the  loop  circulations  would  eventually  converge  as  well  (that  is,  sk  would  tend  towards  zero) 
but  this  usually  required  a  large  number  of  iterations. 

If  a  good  initial  guess  is  available  then  the  number  of  Gauss-Seidel  iterations  is  substantially  reduced.  A 
good  initial  guess  for  attached  flow  can  come  from  another  angle  of  attack.  For  separated  flow,  the  set  of 
loop  circulations  from  the  previous  wake  iteration  provides  a  good  initial  guess  as  well.  Frequently  the 
number  of  iterations  was  large  enough  to  make  Gauss-elimination  the  preferred  scheme. 
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Control  point  location  influences  the  structure  of  |A)  and  the  resulting  loads  as  well.  Use  of  average 
control  points  does  not  necessarily  result  in  [A]  being  strictly  diagonally  dominant.  Nonetheless,  the  diagonal 
elements  are  always  greater  in  absolute  value  than  any  of  the  off-diagonal  elements.  For  the  long  thin  panels 
near  the  nose  of  the  body,  however,  the  diagonal  elements  can  overwhelm  the  off-diagonal  elements,  which 
seems  to  make  the  matrix  product  |A]IG]  sensitive  to  small  changes  in  |G|.  Using  optimal  control  points 
minimizes  the  diagonal  elements  in  |A|  while  leaving  the  off-diagonal  elements  essentially  the  same. 
Numerical  examples  in  Section  4.0  illustrate  the  effect  of  control  placement  on  the  calculated  loads. 

Another  source  of  ill-conditioning  is  abrupt  changes  in  the  lattice  spacing.  A  nearly  regularly  spaced 
uniform  lattice  with  gradual  changes  in  panel  size  and  aspect  ratio  is  recommended  for  reliable  and  efficient 
solutions. 

Computer  resources  also  impose  limitations.  If  an  iterative  scheme  is  used,  the  coefficient  matrix  can  be 
stored  on  a  low -speed  deuce  (disk  or  tape)  and  blocks  of  rows  retrieved  when  needed.  If  a  direct  method  is 
used,  however,  the  coefficient  matrix  is  modified  during  the  solution  process.  If  the  coefficient  matrix  cannot 
be  stored  in  memory .  then  successive  rewrites  of  the  matrix  to  a  low-speed  device  are  necessary.  This  is 
time  consuming  for  a  large  system  of  equations  and  if  pivoting  is  used,  the  penalties  of  having  to  retrieve 
the  coefficient  matrix  hi  column  and  then  row  order  are  considerable.  If  memory  limitations  are  severe,  then 
an  iterative  method  should  probably  be  used.  However,  many  modem  computer  systems  employ  virtual 
memory  management  and  extremely  large  systems  (4000  equations  or  more)  can  be  stored  in  memory  for  the 
time  necessary  to  obtain  a  solution  of  the  linear  equations.  Then  the  memory  can  be  dynamically  returned  to 
the  operating  system  so  that  it  is  not  dedicated  to  the  program  while  other  calculations  are  performed. 

An  ingenious  scheme  has  been  tested  by  James  and  Clark  (1979)  to  help  reduce  the  effort  required  to 
solve  the  linear  equations  They  exploited  the  predictable  nature  of  an  influence  matrix  by  expressing  the 
rows  and  columns  ol  the  matrix  in  terms  of  either  a  Fourier  or  an  eigenvector  series.  They  then  solved  a 
reduced  set  ol  equations  in  terms  of  the  coefficients  of  the  basis  functions,  enabling  them  to  write  an  approx¬ 
imate  solution  ot  the  complete  system  of  equations.  Results  for  two-dimensional  flow  problems  indicated  that 
considerable  savings  for  three-dimensional  problems  could  possibly  be  achieved.  However,  the  structure  and 
approximation  ot  the  influence  matrix  is  much  more  complicated  in  the  latter  case. 

1  he  choice  ot  solution  technique  tor  the  system  of  equations  depends  on  whether  a  good  initial  guess  is 
available  and  on  the  conditioning  and  the  size  of  the  system.  We  varied  our  choice  of  technique  depending 
on  the  circumstances 
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Once  the  no  penetration  condition  on  the  body  is  satisfied,  we  position  the  wake  segments  so  that  they 
are  force -free.  Since  the  influence  ot  each  wake  segment  is  felt  everywhere,  mov  ing  any  of  the  wake 
segments  changes  the  force  tree  position  of  all  the  other  wake  segments  and  results  in  the  no-penetration 
condition  being  violated  on  the  body  While  this  influence  is  minimal  for  those  segments  faraway,  the 
influence  ol  segments  nearby  is  substantial  Thus,  the  order  of  adjustment  of  the  segments  is  important.  We 
tound  that  positioning  the  segments  m  turn  resulted  m  faster  convergence  than  positioning  them  all 
simultaneously  lor  a  fixed  set  of  loop  circulations  Typically,  all  the  segments  in  the  wake  must  be  reposi¬ 
tioned  three  or  lour  times  before  a  reasonably  force  tree  wake  is  obtained.  In  practice,  the  loads  on  the  body 
are  siirpiisingly  msensiti'  e  to  the  wake  position  once  the  wake  has  begun  to  roll  up. 


The  physics  of  the  problem  suggests  judicious  choices  in  the  order  of  wake  segment  adjustment.  If  the 
flow  is  symmetric,  we  adjust  corresponding  pairs  of  segments  on  either  side  of  the  plane  of  symmetry 
simultaneously.  Our  experience  has  been  that  it  is  best  to  adjust  the  furthest  upstream  segments  first,  because 
their  position  seems  to  converge  rather  quickly.  We  do  not  adjust  the  position  of  the  segments  in  the  last  two 
or  three  wake  lines.  Their  converged  position  is  difficult  to  find  and  the  effect  on  the  region  of  interest  is 
minimal.  Evidently,  the  method  tends  to  overestimate  the  distance  the  segments  should  be  moved  from  one 
iterate  to  the  next  since  relaxing  the  predicted  changes  in  the  positions  tends  to  stabilize  the  calculations  and 
usually  results  in  fewer  iterations.  In  fact,  relaxation  factors  as  low  as  one-tenth  are  not  unreasonable, 
especially  in  the  early  stages  of  the  iteration  when  few  numerical  quantities  are  close  to  their  converged 
values. 

Our  model  of  the  wake  corresponds  to  tracking  the  path  of  vorticity  downstream  with  an  Euler  time  in¬ 
tegration  scheme.  Zorea  and  Rom  (1978)  used  both  an  Euler  method  and  a  second-order  Runge-Kutta 
method  to  solve  for  the  wake  position,  but  they  did  not  comment  on  the  differences  they  may  have  en¬ 
countered  between  the  two  methods.  Such  a  scheme  would  certainly  take  more  computer  time  but  would  also 
be  more  stable  than  using  the  lower-order  Euler  method. 

Finding  a  converged  wake  position  with  vortex-lattice  methods  sometimes  requires  several  trial  iterative 
schemes  in  combination  with  under-relaxing  the  wake  motion.  In  fact,  in  solving  for  the  separated  flow  over 
thin  low  aspect-ratio  delta  wings  at  moderate  angles  of  attack,  Mook*  observed  that  certain  combinations  of 
lattice  parameters  could  induce  a  cyclic  behavior  of  wake  position  iterates  after  many  iterations  while  the 
loop  circulations  had  already  converged.  Moreover,  averaging  successive  wake  iterates  only  caused  a  new  set 
of  cyclic  wake  iterates  to  develop.  Under-relaxing  the  motion  of  the  wake  segments  during  the  iteration  pro¬ 
cess  alleviated  the  problem. 

Aparinov,  Belotserkovskii,  Nisht,  and  Solklova  (1976)  claimed  that  the  failure  of  vortex-lattice  methods 
to  converge  on  a  wake  shape  indicated  vortex  breakdown.  This  is  improbable  since  so  many  parameters  in 
the  method  influence  the  sequence  of  iterates.  Moreover,  for  some  problems,  one  combination  of  iteration 
parameters  may  induce  convergence  while  another,  equally  plausible  combination,  may  not. 

3.8  CALCULATION  OF  LOADS 

We  use  the  steady  Bemoullii  equation  to  calculate  the  pressure  coefficient  on  the  body: 


where  Poo 
P 


e 

v 


is  the  freestream  dynamic  pressure 
is  the  local  pressure 
is  the  undisturbed  freestream  velocity 
is  the  fluid  density 
is  the  surface  velocity 


*As  observed  by 


Mook  m  personal  communication 
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We  calculate  the  pressure  coefficient  on  the  body  at  the  control  points.  The  local  normal-force  coeffi¬ 
cient  for  a  fixed  axial  location  on  the  body  is  defined  as: 


C 


n 


cos  0  d0 


(3.10) 


where  r  is  the  radius  of  the  body  and  0  is  the  circumferential  angle  measured  from  the  windward  meridian. 

The  total  velocity  on  the  bound-vortex  sheet  is  the  sum  of  the  contributions  due  to  the  freestream,  all  of 
the  vortex  segments  in  the  (low,  and  one-half  the  jump  in  velocity  across  the  vortex  sheet.  We  calculated 
this  jump  as: 


❖  = 


x  w 


(3.11) 


where  n  is  the  unit  nomial  to  the  sheet  and  w  is  the  vorticity  vector  for  the  panel  as  defined  in  Equation  (3.1). 
3.9  SUMMARY 

In  this  section,  we  have  described  the  vortex-lattice  model  that  replaces  the  continuous  vortex-sheet  model 
presented  in  Section  2.0  and  an  iteration  procedure.  In  the  next  section,  we  present  numerical  examples  that  il¬ 
lustrate  many  of  the  concepts  introduced  here. 


4.0  NUMERICAL  EXAMPLES 


To  demonstrate  the  vortex-lattice  method,  we  treat  the  flow  over  an  inclined  ogive-cylinder  with  a  nose 
length  of  three  diameters  The  circumferential  angle  and  axial  starting  location  of  the  separation  line  are 
assumed  to  the  known 

In  cases  with  attached  flow  we  compare  our  results  with  those  of  the  XYZPF  computer  code  described 
in  Dawson  and  Dean  (1972).  The  XYZPF  code  employs  Hess  and  Smith’s  (1962)  source-distribution  method 
to  calculate  flows  over  nonlifting  bodies  It  is  important  to  note  that  the  same  panel  arrangement  has  been 
used  in  the  calculations  for  both  the  vortex-lattice  method  (VLM)  and  the  source-distribution  method  (SDM). 
All  results  of  the  VLM  and  the  SDM  citied  are  from  calculations  specifically  performed  for  use  in  this 
report 

In  case  with  separated  How  we  compare  our  results  with  those  of  the  concentrated  vortex  core  and 
feeding  sheet  model  of  Hr, son  (1959)  as  reported  by  Jepps  (1977).  Jepps  treated  the  separation  line  on  the 
body  as  the  locus  ot  crossflow  plane  stagnation  points.  He  also  fixed  the  circumferential  separation  angle  as 
we  do  here 

We  measure  the  circumferential  ingle  from  the  windward  meridian.  The  circumferential  angle  of  the 
separation  line  is  called  i  and  the  angie  of  attack  is  a.  All  of  the  cases  considered  here  are  symmetric,  so 
results  are  shown  for  one  half  of  the  body  only.  A  lattice  with  n  panels  on  the  nose  in  the  axial  direction 
and  m  panels  everywhere  on  the  body  in  the  circumferential  direction  is  specified  as  an  “n  x  m"  lattice. 
These  quantities  are  illustrated  in  Figure  4-1. 


Figure  4-1  Definition  of  Lattice  Quantities 

If  necessary,  a  cosine  circumferential  spacing  ot  the  panels  is  used  and  is  adjusted  so  that  the  separation 
line  location  I. ills  on  a  longitudinal  line  in  the  bound-vortex  lattice.  That  is,  the  spacing  in  the  circumferen¬ 
tial  direction  is  related  to  (he  difference  in  the  cosines  of  equal  increments  ol  an  angle.  This  provides  a 
smooth  variation  in  the  size  ol  ad|aeent  panels  If  possible,  a  uniform  circumferential  distribution  of  panels  is 
used  The  circumferential  spacing  is  arranged  to  be  symmetric  about  the  midplane  of  the  body 


For  comparisons  of  results  from  different  calculations  and  numerical  with  experimental  results,  quantities 
are  rarely  available  at  the  common  stations  on  the  body.  To  circumvent  this  problem,  we  use  bicubic  splines 
to  interpolate  the  results.  The  trapezoidal  rule  is  used  for  integrating  the  pressures  around  the  body  to  obtain 
the  normal-force  coefficient. 

4.1  ATTACHED  FLOW  RESULTS 

We  first  examine  the  effect  of  using  area-weighted  branch  circulations  in  the  pressure  calculation  when 
the  lattice  has  an  irregularity  in  the  axial  spacing  in  Figure  4-2  the  planviews  of  two  lattices  are  shown. 

OGIVE  CYLINDER  JUNCTION 
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Figure  4-2  Planview  of  Uniform  and  Nonuniform  Lattices 

Lattice  a  has  uniform  axial  spacing  over  the  entire  length  of  the  body,  while  lattice  b  has  a  sudden  change  in 
the  axial  spacing  immediately  alt  ol  the  ogive-cylinder  function.  The  lattices  have  the  same  number  of 
elements  in  the  circumferential  direction  and  they  are  evenly  distributed. 

In  Figure  4  4  the  pressure  coefficients  calculated  for  the  case  of  zero  angle  of  attack  for  both  lattices  are 
plotted  with  the  expermienial  results  cited  in  Faulkner.  Hess,  and  Giesing  (1964).  Optimal  control  points 
were  used  in  all  the  VI  M  results  shown  The  results  obtained  using  equally-weighted  circulations  near  the 
nose  lor  the  unilorm  and  noimmtorm  lattices  were  virtually  identical  and  are  not  shown.  As  demonstrated  by 
the  figure,  the  use  of  area  weighted  circulations  significantly  reduced  the  evidence  of  the' lattice  irregularity 
in  the  pressure  distribution  I  he  panels  that  join  to  form  the  nose  of  the  body  are  triangular  Therefore,  we 
expect  to  see  differences  in  the  results  using  the  two  calculation  rules  near  the  nose,  as  the  figure 
demonstrates.  The  use  of  equally  weighted  circulations  near  the  nose  resulted  in  pressure  coefficients  closer 
to  experimental  values  Thus,  the  use  of  area  weighted  circulations  provided  less  accurate  results  near  the 
nose  hut  more  accurate  results  near  the  lattice  irregularity  F  his  probably  indicates  that  the  pressure  calcula¬ 
tion  for  panels  that  border  other  panels  of  widely  different  shape  requires  special  consideration.  That  is.  the 
assumptions 
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regarding  vorticity  modelling  used  in  deriving  the  pressure  calculation  scheme  we  employ  are  not  valid  when 
panels  of  widely  different  shapes  are  adjacent.  We  see  virtually  no  difference  among  any  of  the  variations 
used  here  from  just  aft  of  the  nose  to  just  forward  of  the  junction  because  the  panels  do  not  vary  drastically 
from  one  to  the  next. 


Figure  4-3  --  Effect  of  Area-Weighting  Branch  Circulations  with  an  Irregularity  in  the 
Bound- Vortex  Lattice  with  Experiments  of  Faulkner,  Hess,  and  Giesing  (1964) 

The  figures  that  follow  compare  the  results  obtained  by  the  SDM,  the  VLM,  and  experimental  evidence. 
All  calculations  were  performed  on  the  same  computer  with  precisely  the  same  panel  arrangement,  where 
indicated. 

Figure  4-4a  compares  the  results  of  the  SDM,  the  VLM  using  optimal  control  points,  and  experiments. 
As  expected,  the  results  calculated  with  the  SDM  method  agree  more  favorably  with  the  experimental  data 
near  the  nose  of  the  body  than  do  the  results  of  the  VLM.  The  results  of  the  two  methods  are  in¬ 
distinguishable  on  the  cylindrical  portion  of  the  body. 

The  influence  of  the  axial  density  of  elements  on  the  results  of  the  VLM  is  shown  in  Figure  4-4b  for  the 
case  of  zero  angle  of  attack.  An  increase  in  the  axial  density  of  the  elements  provided  for  better  agreement 
with  experimental  data  in  the  nose  region.  A  similar  comparison  for  the  SDM  showed  no  discrenible  dif¬ 
ference.  The  results  suggest  that  the  VLM  is  somewhat  more  sensitive  to  axial  density  of  elements  than  the 
SDM  is  However,  this  sensitivity  is  not  extreme. 

The  influence  of  circumferential  panel  density  on  the  pressure  coefficient  for  the  case  of  zero  angle  of 
attack  is  shown  in  Figure  4  5a  for  the  VLM  and  in  Figure  4-5b  for  the  SDM.  Surprisingly,  the  results  of  the 
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AXIAL  DISTANCE  (diameters) 


f  igure  4-4  Pressure  Coefficient  versus  Axial  Distance  for  a  =  0°.  (a)  Calculated  with  the  VL\1 
and  the  SDM;  and  (h)  the  Effect  of  Axial  Panel  Density  on  the  VLM.  Both  are 
Compared  with  Experimental  Results  of  Faulkner,  Hess,  and  Giesing  (1964) 
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f  igure  4  5  fitted  of  C  ircumferential  Density  of  Elements  (a)  in  the  VLM.  and  (b) 
in  the  SDM.  for  a  -  0°.  as  Compared  with  Experimental  Results  of  Faulkner, 

Hess,  and  Giesing  (1964) 
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VLM  using  4  circumferential  panels  are  more  accurate  than  those  using  14  circumferential  panels.  This  is 
because  in  the  latter  case  the  panels  have  become  long  and  thin.  The  closer  the  panels  are  to  being  square 
the  better  the  results  tend  to  be  due  to  the  assumptions  made  in  the  development  of  the  surface  voticity 
model.  Using  area-weighted  circulations  will  not  alleviate  the  problem  because  adjacent  panels  have  approx¬ 
imately  the  same  area.  The  SDM  seems  to  show  the  reverse  behavior.  However,  the  influence  of  cir¬ 
cumferential  panel  density  on  the  calculated  results  is  markedly  different  when  the  body  is  inclined  to  the 
flow  as  we  illustrate  below. 

The  data  suggest  that  the  results  of  the  VLM  using  optimal  control  points  and  the  results  of  the  SDM 
are  relatively  insensitive  to  panel  density  for  the  case  of  zero  angle  of  attack  so  long  as  the  panels  remain 
approximately  square. 

We  now  investiagte  the  effect  of  control-point  location  on  the  results  of  the  VLM  for  the  case  of  zero 
angle  of  attack.  In  Figure  4-6,  the  pressure  coefficient  calculated  by  the  VLM  using  the  two  choices  of 


Figure  4-6  —  Effect  of  Control  Point  Location  in  the  VLM  for  a  =  0°.  as  Compared 
with  Experiments  of  Faulkner,  Hess,  and  Giesing  (1964) 

control-point  placement  is  shown.  Here,  the  use  of  optimal  control  points  significantly  improved  the  accuracy 
of  the  calculated  pressure  near  the  nose.  The  influence  of  control-point  location  diminishes  as  we  move 
back  along  the  body  because  the  location  of  the  two  control  point  choices  grow  closer  together.  The  two 
control  point  locations  are  coincident  on  the  cylindrical  portion  of  the  body  since  the  panels  there  are  rec¬ 
tangular  in  shape. 


However,  influence  of  control  point  placement  is  entirely  different  when  he  body  is  inclined  ti  the 
flow.  Moreover,  the  influence  varies  with  the  axial  station  and  circumferential  panel  uchsu;.  i  .6ure  4-7 
shows  the  results  for  the  first  four  axial  calculation  stations  for  the  case  of  15°  angle  of  attack  using  the 
VLM  with  optimal  control-point  placement.  The  influence  of  circumferential  panel  density  is  more  pro¬ 
nounced  than  it  was  in  the  case  of  zero  angle  of  attack.  This  influence  diminished  as  we  move  back  along 
the  body. 
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Station  3  xld  =  0.5631 
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Station  2  x/d  =  0.3730 


a  =  15°.  Opt.  CP  s 

-  VLM.  14  x  6 

- VLM.  14  x  10 

- VLM.  14  x  14 


Station  4  xld  =  0.7625 


Figure  4-7  —  Effect  of  Circumferential  Density  of  Elements  in  the  VLM  "using 
Optimal  Control  Points  for  a  =  15° 

In  Figure  4-8.  we  show  the  corresponding  results  using  average  control-point  placement.  The  sensitivity 
to  circumferential  panel  density  is  much  less  severe.  The  reason  that  the  station  axial  locations  are  different 
in  the  two  figures  is  that  average  and  optimal  control  points  are  located  at  different  axial  stations.  Rather 
than  introduce  interpolation  errors  that  might  contaminate  the  comparison,  we  chose  to  present  data  at  the 
stations  used  in  the  calculations. 
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Station  1  x/d  =  0  14 


Station  2  x/d  =  0.32 


Station  3  x/d  =  0.54 


a  =  15°.  Avg  CP's 

-  VLM,  14  x  6 

- VLM.  14  x  10 

- VLM.  14  x  14 

Station  4  x/d  =  0.75 


Figure  4  8  —  lifted  of  Circumferential  Density  of  Elements  in  the  VLM  using 
Average  Control  Points  to.  a  =  15° 

The  corresponding  figure  for  the  SDM  is  presented  in  Figure  4-9.  The  SDM  shows  very  little  sensitivity 
to  circumferential  panel  density  for  this  case. 

Figure  4-10  shows  representative  results  for  all  three  calculations  interpolated  to  common  axial  stations 
for  comparison.  The  difference  between  the  results  of  the  SDM  the  and  VLM  using  optimal  control  points 
grow  smaller  as  we  move  axially  along  the  body,  whereas  the  results  of  the  VLM  using  average  and  optimal 
control  points  grow  more  parallel.  This  is  illustrate-.'  further  in  Figures  4-11  and  4-12. 

The  results  of  the  calculations  are  compared  with  experimental  values  reported  by  Tinling  and  Allen 
(|9f>2)  in  Figure  4-11.  The  experimental  data  were  taken  at  Reynolds  numbers  of  440.000  and  3,000.000 
based  on  the  test  model's  length  of  10  diameters.  At  x/d  =  0.5,  the  pressure  distribution  calculated  with  the 
vortex-lattice  method  using  optimal  Control  points  agrees  better  with  the  experiments  than  the  one  calculated 
using  average  control  points.  The  results  of  calculations  using  average  control  points  appear  almost  parallel 


30 


NORMAL  FORCE  COEFFICIENT 


0.0  1 - 1 -  t-— - J - 1 

0  2  4  6  8 


AXIAL  DISTANCE  (diameter*) 


Figure  4-24  -  Calculated  Attached  and  Separated  Normal-Force  Coefficients  with  Calculations 
of  Jepps  ( 1977)  and  experimental  Results  of  Tinling  and  Allen  (1962)  for 
a  -10".  15".  and  20  for  ^  =  149. 1" 


Figures  4-23  and  4-24  display  tht-  normal-force  coefficient  versus  axial  distance  for  different  separation 
line  locations  and  angles  of  attack.  Both  the  cross-flow  analogy  and  vortex-lattice  method  calculated  forces 
that  level  off  on  the  afterbody  quickly  at  the  10°  angle  of  attack  but  not  at  all  for  a  20°  angle  of  attack.  The 
sharp  peaks  in  the  results  for  the  vortex-lattice  method  are  present  because  of  the  sudden  departure  in  the 


AXIAL  DISTANCE  (diamatera) 

Figure  4-23  —  Calculated  Attached  and  Separated  Normal-Force  Coefficients  with  Calculations  of  Jepps 
(1977)  and  Experimental  Results  of  Tinling  and  Allen  (1962)  for  a  =  10°.  15°,  and  20”  for  =  140° 
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Figure  4-22  —  Calculated  Attached  and  Separated  Pressure  Coefficient  with  Calculations 
of  Jepps  (1977)  and  Experimental  Results  of  Tinling  and  Allen  (1962)  for 
x/d  =  4.5.  6.0.  and  7.5  for  a  =  20°;  +  =  149.1° 


Figure  4-20  Calculated  Attached  and  Separated  Pressure  Coefficient  with  Calculations 
of  Jepps  ( 1077)  and  Hxpenmental  Results  of  Tinting  and  Allen  (1962)  for 
v'd  4.5.  6.0.  and  7.5  for  a  =  10°;  $  =  149.1° 
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Figure  4-19  —  Calculated  Attached  and  Separated  Pressure  Coefficient  with  Calculations 
of  Jepps  (1977)  and  Experimental  Results  of  Tinting  and  Allen  (1962)  for 
x/d  ~  4.5.  6.0,  and  7.5  for  a  -  20”;  41  =  190° 
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Figure  4-18  —  Calculated  Attached  and  Separated  Pressure  Coefficient  with  Calculations 
of  Jepps  (1977)  and  Experimental  Results  of  Tinling  and  Allen  (1962)  for 
x/d  =  4.5,  6.0,  and  7.5  for  a  =  15°;  <{»  =  140° 
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Figure  4-15  —  Calculated  Attached  and  Separated  Flow  Pressure  Coefficient  for 
Calculation  Stations  9  through  12  for  a  =  15° 
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Figure  4-16  —  Calculated  Attached  and  Separated  Flow  Pressure  Coefficient  for 
Calculation  Stations  13  through  16  for  a  =  15° 


Figure  4-14  shows  the  front  view  of  the  same  wakes.  In  the  Figure  all  the  wake  lines  have  been  collaps¬ 
ed  onto  the  viewing  plane.  The  diagram  does  not  represent  a  section  view  of  the  wake  but  rather  all  the 
wake  lines  are  shown  together  as  they  would  be  seen  by  an  observer  in  front  of  the  body  looking  aft.  While 
most  of  the  lines  roll  up  in  a  more  or  less  organized  vortex  core,  a  few  of  the  lines  extend  beyond  the  core 
envelope.  These  lines  originate  near  the  end  of  the  body.  Finding  their  converged  position  proved  to  be  very 
time  consuming  while  at  the  same  time  their  influence  on  the  loads  in  the  region  of  interest  on  the  body  was 
minimal.  For  the  sake  of  economy  in  the  examples  shown,  the  last  three  wake  lines  were  not  adjusted  to  a 
force-free  position,  and  thus  the  few  lines  preceding  them  converged  to  unrealistic  shapes.  The  influence  on 
the  loads  on  the  body  in  the  region  of  interest  was  still  negligible. 


ci  ■  io*  a  •  15”  a  •  20’ 

Figure  4-14  -  Front  View  of  Calculated  Vortex  Wakes  for  a  =  10°.  15°.  and  20° 

To  demonstrate  the  behavior  of  the  circumferential  pressure  distribution  at  different  axial  stations  in  the 
attached  to  the  separated  flow  region  we  compare  the  calculated  attached  and  separated  circumferential 
pressure  distributions  for  several  consecutive  axial  calculation  stations  in  Figures  4-15  and  4-16.  The  Figure 
illustrates  circumferential  and  longitudinal  changes  in  the  pressure  distribution  as  we  move  from  attached  to 
separated  regions  on  the  body. 

The  First  vortex  line  is  attached  to  the  body  between  stations  12  and  13.  Only  a  few  stations  later,  the  shape 
of  the  circumferential  pressure  distribution  does  not  change  from  one  station  to  the  next.  The  pressure  of  the 
separated  vortex  sheet  increases  the  pressure  in  the  front  of  the  separation  line  over  the  value  it  had  for  attached 
(low.  while  the  pressure  behind  the  separation  line  falls  below  the  value  for  attached  flow.  Since  there  arc  only 
two  circumferential  panels  behind  the  separation  line,  the  calculations  do  not  repnxluce  flow  details  there. 

Figures  4-17  through  4-22  display  the  various  circumferential  pressure  distributions  for  10°,  15°,  and 
20"  angles  of  attack  for  x  d  =  4.5.  6.0,  and  7.5  In  all  cases  the  cross-flow  analogy  method  of  Jepps  over¬ 
predicted  the  rise  in  the  pressure  in  front  of  the  separation  line  and  the  results  show  a  sharper  drop  in  the 
pressure  behind  the  separation  line  than  do  those  of  the  present  method.  Neither  method  does  well  behind  the 
separation  line  Both  methods'  results  are  least  accurate  at  the  higher  angles  of  attack.  Note  that  the 
Reynolds  number  effect  on  the  experimental  data  is  more  pronounced  at  higher  angles  of  attack. 
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both  to  those  using  optimal  control  points  and  to  the  experimental  results.  The  difference  in  the  experimental 
results  due  to  Reynolds  number  effects  is  greatest  on  the  leeward  side  of  the  body— precisely  where  the 
results  of  the  VLM  also  show  the  vary  furthest  from  parallel.  Except  for  the  dip  in  the  average  control-point 
results,  the  difference  between  the  two  VLM  calculations  at  this  station  represents  approximately  a  constant 
shift  in  the  pressure  coefficient  throughout  the  circumference  of  the  section.  Such  a  shift  in  the  pressure  coef¬ 
ficient  will  not  influence  the  resultant  normal-force  coefficient  since  the  integral  of  a  constant  over  the  cir¬ 
cumference  of  a  section  vanishes.  Subsequent  boundary-layer  calculations  would  remain  largely  unaffected 
since  the  pressure  gradient  is  also  unaffected  by  a  constant  shift  in  pressure. 

At  x/d  =  2,  (in  Figure  4-1  la)  the  shift  in  the  pressure  distribution  is  smaller  but  still  present.  Also,  the 
results  of  the  VLM  using  optimal  control  points  appear  to  be  slightly  larger  than  the  experimental  results. 

The  results  using  average  control  points  now  even  appear  to  be  slightly  more  accurate  but  this  is  probably 
fortuitous.  There  also  seems  to  be  more  of  a  Reynolds  number  effect  at  this  station. 

In  both  figures,  the  results  of  the  SDM  agree  as  well  as  the  results  of  the  VLM  using  optimal  control 
points. 

Figure  4-12  displays  the  normal-force  coefficient  versus  axial  distance  for  three  different  angles  of  at¬ 
tack.  In  all  three  cases,  results  of  the  VLM  using  average  control  points  agree  better  than  those  of  the  VLM 
using  optimal  control  points,  even  though  the  opposite  behavior  was  experienced  when  viewing  pressure 
distribution  comparisons.  The  reason  for  the  behavior  is  not  obvious  As  expected  the  shift  in  the  pressure 
distribution  due  to  control-point  placement  did  not  substantially  alter  the  normal-force  coefficient  for  the  sta¬ 
tions  we  examined.  Unlortunately  experimental  data  were  unavailable  for  intermediate  stations.  The  results  of 
the  SDM  are  reasonable  for  all  three  cases. 

The  normal-force  coefficient  on  the  afterbody  is  zero  for  attached  flow  because  no  lift  mechanism  is  pre¬ 
sent  The  experimental  data  indicate  that  lift  is  present  aft  of  the  ogive-cylinder  junction  since  the  flow 
begins  to  separate  from  the  body  there. 

The  variation  of  the  experimental  data  with  Reynolds  number  ndicates  that  viscous  effects  are  of  some 
importance.  Including  the  displacement  thickness  of  the  boundary  layer  in  the  dimensions  of  the  body  might 
improve  agreement  between  calculated  and  measured  results 

4.2  SEPARATED  FLOW  RESULTS 

We  now  compare  results  of  the  VLM  in  cases  with  separated  flow  to  the  calculations  of  Jepps  (1977) 
and  the  experimental  results  of  Tinling  and  Allen  (1962).  The  separation  angle  is  adjusted  such  that  the  loca¬ 
tion  of  the  control  point  just  below  (that  is,  just  windward)  of  the  separation  line  is  approximately  the  same 
in  our  calculations  as  it  was  in  Jepps.  Since  we  used  a  circumferential  spacing  that  was  far  coarser  than 
Jepps',  our  separation  angle  is  somewhat  higher  than  Jepps'.  That  is,  the  separation  angle  in  the  present 
calculations  is  located  more  towards  the  leeward  side  of  the  body.  We  chose  Jepps’  work  as  representative  of 
methods  based  on  two-dimensional  analogy. 

In  Figure  4-13,  the  wake  off  one  side  of  the  body  is  shown  for  10°,  15°,  and  20°  angles  of  attack  The 
figure  is  a  computer-generated  drawing  of  the  actual  solutions  calculated  by  the  present  method.  The  separa¬ 
tion  angle  for  all  cases  is  140°  There  are  24  wake  lines  with  the  first  line  having  30  finite  vortex  segments 
and  each  subsequent  line  having  one  segment  fewer  than  the  line  in  front  of  it.  This  allows  for  many  of  the 
wake  lines  to  end  at  approximately  the  same  axial  location. 
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Figure  4-12  —  Calculated  Attached  Normal-Force  Coefficient  Compared  with  the  Experimental 
Results  of  Tinling  and  Allen  (1%2)  for  a  =  10°,  15°,  and  20° 
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lattice  configuration  there  due  to  the  presence  of  the  free-vortex  sheet.  The  bound-vortex  lattice  for  these  ex¬ 
amples  was  tested  to  ensure  that  no  irregularities  in  the  results  were  due  to  irregularities  in  the  bound-vortex 
lattice  spacing.  Jepps  mentioned  also  that  his  method  produced  unreliable  results  in  this  region  of  the  body 
and  he  smoothed  the  normal-force  coefficient  using  values  from  neighboring  stations  to  make  the  results  look 
more  realistic.  We  have  done  no  such  smoothing  here.  Once  again,  Reynolds  number  effects  are  more 
significant  at  the  higher  angles  of  attack. 

The  results  for  the  two  separation  angles  used  are  shown  in  Figure  4-25  for  the  case  of  a  15°  angle  of 
attack.  It  is  interesting  that  the  change  in  separation  line  location  caused  an  approximately  constant  shift  in 
the  normal-force  coefficient  on  the  lifting  portion  of  the  body.  A  shift  in  the  separation  line  circumferential 
angle,  then,  has  an  effect  on  the  loads  that  is  similar  to  the  effect  of  a  change  in  the  Reynolds  number. 


Figure  4-25  —  Calculated  Separated  Normal-Force  Coefficient  for  a  =  15°,  ^  =  140°  and 
149.1°  with  Experimental  Results  of  Tinling  and  Allen  (1962) 
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Figure  4-26  shows  the  normal-force  coefficient  versus  axial  distance  along  the  body  for  each  iteration  of 
the  loop  circulations.  After  only  three  iterations  the  normal-force  coefficient  has  almost  converged. 

We  plotted  the  loop  circulation  strength  versus  circumferential  angle  for  stations  5  and  25  in  Figure 
4-27.  Station  5  is  located  on  the  nose  of  the  body  and  the  flow  is  fully  attached.  Station  25  is  located  on  the 
afterbody  and  separation  is  present  there.  The  shape  of  the  curves  is  similar  in  both  cases  in  the  attached 
region  of  the  body.  That  is,  the  windward  side  of  the  body  experiences  a  similar  change  in  the  loop  circula¬ 
tion  strength  as  we  move  circumferentially  around  the  body.  Since  the  branch  circulations  (and  hence  the 
velocity  field)  are  dependent  on  the  shape  of  the  curve,  we  could  fix  the  loop  circulations  in  the  attached 
flow  regions  of  the  body  and  calculate  only  those  loop  circulations  which  are  significantly  affected  by  the 
presence  of  the  separated  wakes.  This  would  dramatically  reduce  the  computer  time  needed  for  a  solution.  It 
would  cause  an  error  in  the  branch  circulations  of  the  segments  which  border  an  opening  on  the  body  sur¬ 
face,  but  typically  the  only  opening  is  perhaps  an  open  tail,  which  would  have  a  minimal  effect  on  the 
upstream  portion  of  the  flow  field. 


Figure  4-26  —  Convergence  of  Normal-Force  Coefficient  with  Wake  Iteration 
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5.0  CONCLUSIONS  AND  RECOMMENDATIONS 


We  developed  a  nonlinear  vortex-lattice  method  to  treat  the  steady  flow  over  prolate  bodies  with  open 
separation  moving  through  an  inviscid  incompressible  fluid.  The  strength  and  position  of  the  body  wake  is 
found  as  part  of  the  solution.  Specifically,  flows  with  smooth-surface  separation  have  been  considered  in 
contrast  to  the  flows  with  sharp-edged  separation  that  have  been  treated  with  vortex-lattice  methods  in  the 
past . 

From  our  comparisons,  we  can  conclude  that  the  vortex -lattice  method  is  more  sensitive  to  lattice  ar¬ 
rangement  than  the  source-distribution  method  of  Hess  and  Smith.  We  found  that  use  of  area-weighted  cir¬ 
culations  significantly  reduced  the  effect  of  lattice  irregularities  in  the  calculated  results.  Use  of  optimal  con¬ 
trol  points  with  the  vortex-lattice  method  led  to  improved  results  for  the  zero  angle  of  attack  case  and  worse 
results  when  the  body  was  at  incidence.  Moreover,  using  optimal  control  points  made  the  results  more  sen¬ 
sitive  to  lattice  design.  Also,  we  found  that  if  the  panels  become  long  and  thin  in  one  direction,  the  accuracy 
of  the  results  can  be  degraded  even  though  a  denser  lattice  arrangement  is  being  used. 

In  the  case  of  separated  flow,  we  found  that  both  the  vortex-lattice  method  and  the  two-dimensional 
cross-analogy  Bryson  (1959)  as  applied  by  Jepps  (1977)  overpredict  the  pressure  in  front  of  a  separation  line 
and  underpredict  the  pressure  behind  a  separation  line.  Neither  method  predicts  accurate  loads  near  the 
beginning  of  a  separation  line. 

While  we  considered  only  steady  symmetric  flows,  vortex-lattice  methods  are  applicable  to  asymmetric 
unsteady  flows  as  well.  However,  an  unsteady  separation  line  would  be  necessary.  The  motion  of  a  separa¬ 
tion  line  on  a  body  during  an  unsteady  maneuver  may  be  difficult  to  model. 

The  general  trend  of  the  results  strongly  indicate  that  the  accurate  use  of  the  vortex-lattice  method  re¬ 
quires  great  care  in  the  choice  of  lattice  parameters-particularly  in  the  areas  of  panel  arrangement  and  con¬ 
trol  point  placement. 

For  immediate  areas  of  study  we  recommend  the  following: 

1  Modification  of  the  lattice  to  allow  the  free-vortex  sheet  to  leave  the  body  tangentially. 

2.  Exploration  of  means  to  reduce  the  deviations  in  the  pressures  near  the  beginning  of  a  free-vortex 
sheet  and  behind  it. 

3.  Investigation  of  paneling  techniques  that  allow  for  the  separation  line  angle  to  vary  along  the  body. 

4.  Coupling  of  the  method  with  a  boundary -layer  scheme  to  solve  for  the  location  of  the  separation 
lines  on  the  body  via  a  viscous-inviscid  calculation  procedure. 
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